Doping dependence of spin excitations in the stripe phase of high-T c superconductors 
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Based on the time-dependent Gutzwiller approximation for the extended Hubbard model we 
calculate the energy and momentum dependence of spin excitations for striped ground states. Our 
starting point correctly reproduces the observed doping dependence of the incommensurability in 
La-based cuprates and the dispersion of magnetic modes in the insulating parent compound. This 
allows us to make quantitative predictions for the doping evolution of the dispersion of magnetic 
modes in the stripe phase including the energy u>o and intensity of the resonance peak as well as the 
velocity c of the spin- wave like Goldstone mode. In the underdoped regime nu < 1/8 we find a weak 
linear dependence of u>o on doping whereas the resonance energy significantly shifts to higher values 
when the charge concentration in the stripes starts to deviate from half-filling for rij, > 1/8. The 
velocity c is non-monotonous with a minimum at 1/8 in coincidence with a well known anomaly in 
T c . Our calculations are in good agreement with available experimental data. We also compare our 
results with analogous computations based on linear spin-wave theory. 

PACS numbers: 71.28.+d,71.10.-w,74.72.-h,74.25.Ha 



I. INTRODUCTION 

The development of a theory for high-T c superconduc- 
tivity presupposes a detailed understanding of the normal 
state collective excitations in the charge and magnetic 
channels. This paper is devoted to the latter channel 
where we explore the spectrum of spin fluctuations based 
on a stripe scenario extending the results presented in a 
previous short communication. 1 

The large variety of cuprate materials has often nar- 
rowed the view on the properties of individual com- 
pounds rather than the common features between them. 
In YBa2Cu30j, (YBCO) compounds inelastic neutron 
scattering (INS) experiments have usually emphasized 
the "resonance mode"^* 3 ^*^, an inelastic feature that 
sharpens as a delta-like excitation at frequency loq below 
T c and is located around the commensurate antiferro- 
magnetic wave vector Qaf — (tt, 7t) (we set the lattice 
constant a = 1 and restore it when convenient for clar- 
ity). On the other hand in lanthanum based cuprates 
(LCO) the emphasis has been put on incommensurate 
features, either elastic or low-energy inelastic as for ex- 
ample Goldstone like modes emerging from incommensu- 
rate wave vectorsi 7 ' 8 ' 9 ' 10 : 11 ' 12 : 13 

In LCO the incommensurate low-energy magnetic scat- 
tering occurs at wave vectors Q s = (w ± 2ire, it) and 
Q s = (tt, 7r ± 2ire) and e depends linearly on doping 
e = iih up to rih ~ l/8m& Here rih is the number of 
added holes per planar Cu with respect to the parent 
insulating compound. This behavior is compatible with 
stripe like modulations of charge and spin having a linear 
concentration of added holes v = 1/2 (so-called half-filled 
stripes), that is ~ 1 hole per every 2nd elementary unit 
cell along the stripe and with charge (spin) periodicity d 
(2d) perpendicular to the stripe so that e = l/(2d) and 



n h = vj d. 

Regarding superconductivity Wakimoto and 
collaborators^ have recently found a direct relation 
between T c and the intensity of low energy scattering in 
the overdoped regime that together with other studies^ 
establishes a direct relation between superconductivity 
and incommensurate scattering parallel to the Cu-0 
bond across the full phase diagram of LCO compounds 
and makes more stringent the understanding of the 
magnetic excitations. 

Despite the different emphasis mentioned above, some 
investigations in the last years started to focus on com- 
mon features suggesting a universal phenomenology in 
the spin dynamics between different classes of high-T c 
cuprates. These similarities comprise the low-energy in- 
commensurate scattering in YBCOi& which has a similar 
dependence on doping as in LCO, namely a linear in- 
crease at small doping followed by a saturation at larger 
concentrations. Even static incommensurate order has 
been reportecU 7 " in strongly underdoped YBCO. More- 
over, in detwinned samples of YBCO the incommensu- 
rate scattering has been shown to be anisotropic sug- 
gesting underlying stripe correlations^ The presence of 
a "rigid stripe array" has been questioned in the more 
recent Ref. 19 which also reports an anisotropic response 
although in a more subtle way. We will come back to this 
point in Sec. IIVI 

Finally, recent INS measurements show a remarkable 
similarity of the magnetic response of La2- a; Ba 2 ;Cu04 
(LBCO)ii and YBCOSS over a broad range of momenta 
and energy. Both experiments have revealed low en- 
ergy incommensurate excitations which, with increas- 
ing energy, continuously disperse towards the resonance 
mode at Q^iILSE An INS study on optimally doped 
La 2 _ 3 ;Sr 2 ;Cu04 (LSCO)^ 3 - has also resolved the disper- 
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sion of the spin-wave modes emerging from the incom- 
mensurate positions which closely resemble correspond- 
ing data of YBa2Cu306.85>^ 

Above wo the magnetic excitations start to form a ring 
or square shaped pattern around Qaf which increases 
in size with increasing energy. In twinned samples low 
energy acoustic modes together with the high-energy ex- 
citations thus display an 'X'-like dispersion. These fea- 
tures are now well established in both LCO (eventually 
codoped with Nd or Ba)ili4 and YBCO"' 19 ' 2 "' 22 ' 28 ' 24 
while in Bi- 25 and Tl-materials^ only the resonance peak 
has been resolved yet. 

The two most prominent explanations for the features 
in the magnetic scattering are (a) a scenario based on 
dispersing two-particle bound states induced by the an- 
tiferromagnetic (AFM) correlations in a d-wave super- 
conducting syste m 27 i 28 i 29 i 30 i 31 i 32 i 33 i 34 i 35 and (b) theories 
which rely on the presence of stripe correlations in the 
ground state of the systemi 1 ! 36 ! 37 ! 38 ! 39 ! 40 ! 41 ! 42 Approaches 
belonging to (a) are rather popular in describing spin ex- 
citations in the superconducting state of YBCO materi- 
als. The ground state is assumed to be an homogeneous 
metal from the starti 27 i 28 i 29 i 30 i 31 i 32 i 33 i 34 Superconductiv- 
ity on top of this induces a spin-gap on the dynamical 
susceptibility and an RPA resummation of the AFM in- 
teraction creates exciton peaks in the imaginary part of 
the magnetic susceptibility. Due to the special frequency 
and momentum structure of the non-interacting spin re- 
sponse functions in a d-wave superconducting system, the 
strongest weight is acquired by the 'resonance peak' at 
the AFM wave- vector and frequency loq close to the spin- 
gap energy. Above and below ujq also an incommensurate 
response can be obtained the structure of which strongly 
depends on the underlying tight-binding dispersion. This 
'excitonic' picture finds support in the fact that in opti- 
mally doped YBCO the resonance appears below T c and 
offers a natural explanation for the relation loq ~ kT c 
as observed in YBCO (see e.g. Ref. 0). However in un- 
derdoped YBCO the resonance is observed above and 
upon applying this approach to a quantitative description 
of magnetic excitations in LSCO it has been shown to be 
incompatible with the low energy responsei 3 ^ 3 ^ More- 
over, in order to properly account for the doping depen- 
dence of the low energy incommensurability in YBCO 
this approach has to rely on peculiar nesting properties 
of the Fermi surface over a large doping ranged 

The incommensurate correlations observed on several 
compounds have encouraged theories based on the pres- 
ence of stripe correlations in the ground state of the 



system 
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The starting point is a bro- 



ken symmetry state where at least C4 lattice symmetry is 
broken, translation symmetry is usually broken (at least 
in explicit computations) and the breaking of SU(2) spin 
rotational symmetry de pends o n the approach. If SU(2) 
is conserved as in Refs. I36l37ll3^ the system has a spin 
gap in the spectrum whereas if it is broken ^LiML 4 ^ 4 * 4 ^ 
as in the present work, one finds a Goldstone magnetic 
mode emerging from the incommensurate wave vector. 



A spin gap is observed in YBCO compounds whereas a 
Goldstone like mode is observed in LCO. The response 
of both systems is very similar above the gappi 3 - We will 
come back to this point in Sec. [V] 

In our SU(2) broken approach, at small doping, par- 
tially filled stripes are stabl o 46 ' 47 explaining the linear 
dependence of the incommensurability on doping for 
ri/j < 1/8. 48 ' 49 - 50 - 51 Our previous computations in the 
Gutzwiller approximation (GA)£2i£i have shown that due 
to their transverse extension, stripes with d < 4 are ener- 
getically unfavorable so that for rih > 1/8 the incommen- 
surability remains at e ~ 1/8 in accord with experiment 
and the stripe filling starts to increase beyond the value 
v = 1/2. This picture is also in agreement with dynami- 
cal mean-field theory computations! 4 ^* 4 ^ 

Within this scenario the resonance feature appears as 
a saddle-point (or a local maximum) in the dispersion 
relation of magnetic excitations as quantitatively eval- 
uated for LBCO in our previous worki The idea that 
the resonance may be seen as the lowest energy magnetic 
excitation at Qaf in an incommensurate system was al- 
ready proposed in Ref. [52- 

Within the stripe scenario another line of thought has 
appeared based on simplified effective spin only models. 
These studies map the stripe structure to an array of 
coupled n-leg spin ladders with n even. This allows to 
consider states where SU(2) is not broken and which ex- 
hibit a spin gap as in YBCOji 3 ^ 3 ^ 3 ^ 9 . However, for some 
dopings n is expected to be odd and the system should be 
gapless. We are not aware of evidence for this even/odd 
phenomenon. 

In this paper we present detailed computations of the 
dynamical magnetic structure factor on top of the static 
GA stripe textures using the time dependent Gutzwiller 
approximation (TDGA). 53 ' 54 In addition we compare the 
TDGA results with linear spin-wave theory (LSWT) 
computations. We present the doping dependence of 
magnetic excitations extending our previous results 4 - for 
doping rih = 1/8. To be specific we choose parameters 
as appropriate for LBCO or LSCO where our investi- 
gation provides a quantitative prediction for the relation 
between ojo, the spin- wave velocity, the dispersion of high 
energy modes, etc. as a function of doping rih- As in 
our previous report parameters are fixed by fitting the 
dispersion relation in the insulator and by requiring the 
correct filling v of the stripes. For the rest the compu- 
tations in the doped phase can be considered without 
free parameters and are in good agreement with data 
at Uh — l/8i 4 on LBCO. We compare the new results 
with experimental datai 3 - at rih — 0.1 and n/, = 0.16. 
Encouraged by the experimental finding of universality 
of spin excitations in different cuprate materials we also 
attempt to make connection of our results with exper- 
imental features observed in YBCO materials. In this 
case however, due to the presence of the spin gap and 
the bilayer structure, our results will be more qualitative 
than quantitative. 

This paper is organized as follows. In Sec. [H] we present 
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the model and outline the formalism for the computation 
of magnetic excitations. In Sec. Illll we show how the ap- 
propriate parameter set can be obtained from experimen- 
tal and theoretical considerations. In addition we present 
the saddle-point solutions on which the subsequent RPA 
fluctuations are computed. In Sec. HVI we present results 
for the doping dependence of magnetic excitations in the 
stripe phase and discuss these spectra with respect to the 
present experimental situation. Finally we conclude our 
considerations in Sec. |V] 



II. MODEL AND FORMALISM 

Our investigations are based on the one-band Hubbard 
model with hopping restricted to nearest (ty = —t) and 
next nearest (ty = —t') neighbors 



Here c-^ destroys (creates) an electron with spin a at site 
i, and rij jCr = cj ff Cj )er . U is the on-site Hubbard repulsion. 



A. Time-dependent Gutzwiller approximation 

As a starting point we treat the model Eq. JJJ 
within an unrestricted Gutzwiller approximation (GA) 
in the spin-rotational invariant formulationiSi2i*£I Basi- 
cally one constructs a Gutzwiller wave function by 
applying a projector to a Slater determinant \SD) which 
reduces the double occupancy. The Slater determinant is 
allowed to have an inhomogeneous charge and spin dis- 
tribution describing generalized spin and charge density 
waves determined variationallyiS 
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FIG. 1: Unit cell used for the d = 4 stripe computation. The 
length of the spin arrows corresponds to the LSWT compu- 
tations. (For the GA computations the spin magnitude is 
modulated as shown in Fig.0J. We also show the interactions 
used in the LSWT computation: a uniform value of the AFM 
interaction J across all the bonds except for the ferromagnetic 
bonds where the sign is opposite to ensure stability and the 
modulus is Jf = 0.2J. 



The resulting energy functional E GA = 
treated in the Gutzwiller approximation reads as^ 

E ^ ' ijj z i,cri ,<jZj a ( j 2 (c^ gi Cj CT2 )o -f- U S ' Dj, 

i 

(1) 



JJ,(T,CTi,(T2 



where the matrix Zj is given by 



and 



Zi + cos 2 ^ + Zi sin 2 ^ ^r[zi + ~ Zi ]cos$ 
n4r [zi + — Zi~] cos $ Zi + sin 2 $ + zf cos 2 % 



(2) 



tan $ = 



(A 



±\2 



yi-nj + AAf +\?VD~j 
.y(l-A-(A±) 2 )(n,-A-(An 2 )' 
n,/2 - A ± 5fVl + tan 2 $. 



The above expressions are given in terms of the one- 
body density matrix pi a ,jcr' = ( c L c io-')o in the Slater 
determinant \SD). For brevity we denote such averages 
with (...)o and introduce the notations: = (<4»c,-j,)o, 

n * = E<r( c L c ^}o, etc. 

In order to obtain generally inhomogeneous solutions 
(stripes, etc.) one has to minimize Eq. ^ with respect 
to the double occupancy parameters A = (^\ni^nii\^>) 
and with respect to p under the constraint that p corre- 
sponds to a Slater determinant. For simplicity we restrict 
our saddle-point solutions to the limit S i = 0, where the 
matrix is diagonal and the standard Gutzwiller energy 
functional as derived by Gebhard 58 ^ or Kotliar and Ruck- 
enstein is recovered^ (The spin-rotational invariant for- 
mulation is necessary for the fluctuation computation). 

We have previously shown that with realistic param- 
eters this approach leads to vertical, metallic stripes in 
agreement with experimental data^2*£i Here we have re- 
peated these calculations in large systems typically con- 
sisting of N ~ 100 x 100 lattice sites in order to ob- 
tain sufficient momentum resolution for the computa- 
tion of spin excitations. For stripes oriented along the 
y-direction this is most conveniently done by decompos- 
ing the lattice into N ce u s unit cells of size N a = d X 2 
with N = N ce u s N a . The unit cell for the d = 4 solu- 
tion is schematized in Fig. ^ The locations of the unit 
cells define a Bravais "superlattice" iSi The positions of 
the Bravais superlattice are given by R = niai + ni&i 
with rii integer. For even d the Bravais lattice is gen- 
erated by the elementary translations ai = (c?, 1) and 
&2 = (0, 2) whereas for d being odd we have ai = (d, 0) 
and &2 — (0,2). From the translation vectors we can 
obtain the generators of the reciprocal superlattice 
defined by the relation a^.b^ = 27r<5y and construct the 
associated magnetic Brillouin zone (MBZ). For odd d the 
latter is simply a rectangle spanned by the reciprocal su- 
perlattice vectors bi = (27r/d, 0) and \>2 = (0, ir) whereas 
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FIG. 2: (Color online) The central polygon is the first mag- 
netic Brillouin zone for d = 4 stripes. Thick lines enclose the 
extended Brillouin zone. The dots are a set of magnetic recip- 
rocal superlattice vectors that define a set of higher magnetic 
Brillouin zones with total volume equal to the EBZ volume. 
Magnetic modes and single particle states can be classified 
by the momentum in the central magnetic Brillouin zone and 
a band index. The magnetic reciprocal superlattice vectors 
closest to (±7r, ±7r) are referred to as Q s . 



for even d the shape is more complex due to the non- 
orthogonality of the translation vectors. 

As an example we show in Fig. |2 the first magnetic 
Brillouin zone (MBZ) and the extended Brillouin zone 
(EBZ) for d = 4 stripes. The momenta in the EBZ label 
the N plane wave states of the system. Alternatively one 
can label the plane wave states by a reduced momentum 
k in the first MBZ (the central polygon in Fig. |2J) and 
a set of N a reciprocal superlattice vectors Q = m\h\ + 
77i2b2 with rrii integer, the momentum of the plane wave 
state being given by Q + k. The dots in Fig.|2indicate a 
possible choice for this set. Each dot defines also a higher 
MBZ's with the dot at the center. 

Magnetic excitations are obtained by computing 
random-phase approximation (RPA) like fluctuations on 
top of the GA saddle-poinljSSi within the TDGA which 
fulfills standard sum rules and on small clusters yields 
excitation spectra in very good agreement with exact 
diagonalizationiSi 

The basic ingredient of this approach is the interaction 
kernel, fn the spin channel this is obtained by expanding 



Eq. up to second order in the spin fluctuations^ 
SE spm = J^VijSSfSSr 

; Ml[SSt 6(4^)0 + s(cl lCll ) ss. 



(3) 



£ 



+ £ M^SS+Sic]^ + 5(c^ Cjl ) SSr}. 

ij 

Here we have defined the following interaction matrices 

d 2 z, 



V„ 



£ 

3,<r 



%[( c L c j>)o + (c ja c lr7 )o}z[ 



w ^i,a,a 
3 



Vi 



Ml 



dS+dS,, 

est dsr 



t-z- 



BSt 



(4) 

(5) 
(6) 



Basically the TDGA is closely related to the time- 
dependent Hartree-Fock (HF) approximation which in 
the limit of small oscillations corresponds to the RPA. 63 
The structure of the equations is the same but the in- 
teraction kernel is substantially different. Indeed in the 
present approach there is a space dependency of the ef- 
fective on-site term Vu (instead of the bare 'U') which 
follows the modulation of the underlying saddle point so- 
lution (cf. Fig. In addition the GA+RPA formalism 
generates effective interactions between spin fluctuations 
on different sites (~ Vi^j ) and between intersite spin flip 
processes (~ Mg). However, the strength of the latter 
interactions is in general much smaller than the on-site 
contribution Vu. 

In order to evaluate the RPA series we first have to 
evaluate the fzeroth order') susceptibilities on the GA 
level, e.g. 



st,0f~ 
A. a m 



(7) 



- Idte 



f (T^(q + Q„, t)A*(-q - Q m , 0))ga- 



Here the vector {A s } is composed of the following oper- 
ators: A\ = S' + (q+ Q n ) is the momentum representa- 
tion of the local spin fluctuation 5+ = cLc^. A\ +Q = 
T Q (q + Q n ) are the Fourier transformed of the intersite 
spin fluctuations (y^Ci+Ral- Both fluctuations are mixed 
by the term Mg in Eq. © (cf. Ref. |H for details). For 
a model with next- and next-nearest neighbor hopping 
we have a = 1, . . . , 8 orthogonal operators T a (q + Q„) 
so that together with 5 + (q+Q„) the complete zeroth 
order correlation function x«m (Q; w ) is a 18d X 18d ma- 
trix. On the other hand we can also express the Fourier 
transformed interaction Eq. ® in terms of a 18d x 18d 
dimensional matrix T(q) 



§E s P in = £ A s (q + Q n )r^ m (q)A t (-q-Q r . 



(8) 
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so that formally the RPA summation can be written as 
l(q)=^(q)-^(q)E(q)^(q). (9) 

Upon inverting this equation we finally obtain the 
transverse spin correlation function as Xnm(q) = Xnm(q)- 

B. Spin-wave theory calculations 

In order to interpret the GA+RPA results we have 
made a LSWT computa tion of the dispersion similar to 
those performed in Refs. I4ll42l For simplicity we take a 
model with only AFM nearest neighbor magnetic inter- 
actions J equal everywhere except for the ferromagnetic 
bonds where the interaction has opposite sign (to ensure 
stability) and magnitude Jp as shown schematically in 
Fig. HI 

C. Dynamical structure factor 

Neutron scattering experiments probe the longitudinal 
and transverse components of the dynamical structure 
factor S a/3 (q,ui) weighted by polarization factors iSStSl 

S e ff{q,oj) = J2((l-£))do m S aa (q,to) (10) 

a 

= 77|| S 11 (q, u) + V±S ± (q, to) 

where q = k — k' (q = <l/\q\), k (fe') is the initial (final) 
wave vector of neutrons and an average over the orienta- 
tion of domains, (■■■}dom has been done. In Eq. (|10|l we 
have defined the longitudinal (||) and transverse contribu- 
tions (_L) with respect to the ordered magnetic moment. 

Unpolarized experiments are sensitive to relatively 
sharp features which are expected to correspond mainly 
to transverse fluctuations^ therefore in the present study 
we restrict to the latter S e *f « r/xS 2 -. The aver- 
age polarization factor r)± depends on the orientation of 
the magnetic domains. For a distribution of domains 
isotropic in spin-space r/± = 2/3. For a non isotropic 
distribution rj± depends also on the scattering geometry 
with the minimum value r]± = 1/2& 

In terms of the reduced momentum q and the recip- 
rocal superlattice vector Q„ defined in the previous sub- 
section the structure factor is given by, 

S ± (q+Q n) cj) = ^^-Z d Jdt expH^)Im X L(q,t). 

(11) 

Here we have included an intensity renormalization Zd 
similar to LSWT. For a half-filled system with a SDW 
and in the limit of large U, GA+RPA, HF+RPA and 
LSWT give the same results. The latter is well known 
to overestimate the intensities of the spin-wave modes 
roughly by a factor of two and the same problem arises 
in the RPA approaches. GA+RPA, HF+RPA interpo- 
late smoothly between this strong coupling limit and the 



noninteracting case, however, in the latter case Z^ = 1. 
The value appropriate for cuprate parameters will be in- 
termediate between these two limits. We refer to Ref.lo7l 
where Zd and i]± are discussed in detail. 

S (q, to) in general acquires different weight in MBZ's 
centered on different Q„ (cf. Fig. [5J. Since we are not 
considering the magnetic form factor, S ± (q,u>) repeats 
with the periodicity of the EBZ reciprocal lattice vectors 
(0,2tt) and (2tt,0). 

The stripes correspond to a long-range magnetically 
ordered state which within our RPA-like scheme induces 
a Goldstone mode going to zero frequency at momenta 
q = 0. However, this mode has negligible weight at q ~ 
but acquires a significant strength in the higher magnetic 
zones at q + Q with the magnetic reciprocal lattice vec- 
tor Q + 1 0. The largest spectral weight, as shown below, 
is located at the dots closest to (±7r, ±n) in Fig. El These 
points are referred to as the incommensurate wave vector 



III. PARAMETERS AND SADDLE-POINT 
SOLUTIONS 

With respect to our previous investigation 1 we have 
significantly improved the resolution in frequency and 
momentum space which allows us to consider up to 
10.000 lattice sites and frequency intervals of ~ lmeV 
as discussed in the preceding section. Within the present 
one-band description we fix parameters by the following 
two conditions. 

1. ) The value of the nearest neighbor hopping has 
been fixed to t'/t = —0.2 according to Refs. I5lll68l In 
Ref. 0^] we found that the stripe filling is very sensitive 
to t'/t and it turned out that a ratio of t'/t = —0.2 is 
necessary to describe half-filled stripes as inferred from 
the incommensurability. 10 This value is close to the value 
obtained in a local density approximation studypS 

2. ) The spin- wave dispersion of the undoped sys- 
tem along the magnetic zone boundary poses a rather 
stringent constraint on the value of U. In fact, LSWT 
for the Heisenberg model with only nearest-neighbor ex- 
change interactions J = At 2 /U does not yield any dis- 
persion at all. Only higher order in t/U exchange con- 
tributions (cyclic exchange being the most important 
among them£2*22i£i*2£) start to produce a dispersion in 
the magnetic excitations along the magnetic zone bound- 
ary which is therefore very sensitive to the value of U. 

We have computed the magnon spectra for the half- 
filled system within the GA+RPA approach and fitted 
the spin- wave dispersion of undoped La2Cu04 measured 
at T = lOKH As can be seen from Fig. |21the parameter 
set U/t = 7.5 yields excellent agreement with the mea- 
sured dispersion also along the magnetic zone boundary. 
The fitting also fixes the overall energy scale t = 342meV. 
This dispersion is in slightly better agreement with ex- 
periment than a previous one§£ with U/t = 8 which was 
used in Ref. Q] and enhances the already good agreement 
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between theory and experiment for the position of the 
resonance mode at n h = 0.125 discussed below. We note 
that for U/t — 10 the fitting in the insulator is signifi- 
cantly worsens 7 - All our further calculations are without 
free parameters. It should be mentioned that a similar 
estimate for U / t has been obtained from quantum Monte- 
Carlo calculations of the magnon dispersion 73 and is in 
the range U/t ~ 6,8 obtained from mapping the three 
band Hubbard model to the one-band models More- 
over, a similar value of U/t ~ 8 has been obtained in a 
recent fit of the doping dependence of spectral weights in 
optical data^ 

Finally we note that although the present parame- 
ter set was optimized with magnetic properties it also 
gives reasonable results for charge excitations. Within 
GA+RPA one obtains a gap in the optical excitations 
of 1.8 eV which is close to the 2eV value of reflectivity 
experiments^ 

The saddle point solutions have been discussed in 
Ref. |yl Usually we consider stripe textures oriented 
along the y axis (vertical stripes) where site centered (SC) 
and bond centered (BC) solutions have similar energies. 
The magnetic unit cell for BC solutions is shown schemat- 
ically in Fig. ^ Spin are ordered antiferromagnetically 
in the vertical direction whereas the ferromagnetic spin 
alignment between sites at transverse positions and 1 
indicates the presence of a domain wall where also doped 
holes accumulate. In this paper we label stripes by the 
horizontal distance d between the hole rich domain walls, 
the doping and the symmetry (BC or SC). In order 
to compare with experimental data in L^-xSr/zCuC^ or 
La2- 2; Ba a: Cu04 we assume rih = x. 

Within the present one-band model calculation we find 
that the more favorable solutions at small doping are SC 
but with an energy difference to BC stripes that is neg- 
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FIG. 3: Energy and wave vector dependence of magnetic 
excitations in the half-filled system as obtained from the 
GA+RPA approach for U/t = 7.5, t'/t = -0.2 and t = 
342meV. Squares and triangles correspond to data points from 
INS experiments at T — 10 A" on La2Cu04 by Coldea et aljS 
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FIG. 4: Profile of the charge (top panels), staggered magnetic 
order (middle panels) and effective interaction Vu (lowest pan- 
els) for d = 4 and d = 10 bond-centered stripes in the GA. 
The dashed line in the lowest panels indicates the 'bare' spin 



interaction of the HF+RPA VA 
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ligibly small in large systems. This quasi-degeneracy is 
also found in other one-band calculations!^^ 7 ^^ On 
the other hand BC textures turn out to be more sta- 
ble in the more accurate three-band Hubbard models 
and in first principle computations^ Therefore we pay 
special attention to BC solutions at low dopings. In 
the three-band model the quasi-degeneracy appears at 
slightly larger doping and was related to the anomalous 
optical properties^ 

In Fig.^Jwe show a cut through the charge (top panels) 
and staggered magnetic moment profile (central panels) 
for BC stripe solutions and stripe separation d = 4, 10 
evaluated within the GA. For distant stripes (d = 10 
m Fig. EI) the regions between the domain walls resemble 
the spin-density wave solution of the undoped system (i.e. 
n w 1) whereas doped holes progressively populate the 
AFM domains when stripes become closer. The moment 
in the insulating regions is close to 0.8/j.b and is reduced 
with respect to the full value of lfj,B due to covalency. — 
This is larger than the expected value in the insulator 
~ O.GfiB because we are not including the feedback of 
the RPA fluctuations on the single-particle expectation 
values which, for the magnetic moment, is dominated by 
transverse fluctuations^ 

The weight of the magnetic Bragg peaks is determined 
by the Fourier transform of the magnetic moment in the 
unit cellAS 7 . From the more intense Bragg peak one can 
define an "ordered moment" . For LCO codoped with 
Nd it was found£ that this is about O.I^b whereas the 
present static GA approach yields ~ 0.4/is. 67 Apart from 
the lack of transverse fluctuations which may be even 
more relevant in the stripe phase it is likely that disorder, 
not taken into account in our investigation, will depress 
the weight of Bragg peaks since they require coherence 
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across many unit cells. In this regard it is interesting 
that estimates from local probes, which are less sensi- 
tive to disorder, are much larger and amount to ordered 
moments of ~ This should be compared with 

the average of the absolute value of the static moment 
m GA 1=3 0.56/i,b obtained within the GA. In this case the 
overshoot is expected due to the omission of transverse 
fluctuations. 

Obviously the quality of our approach relies on the 
presence of stripe correlations in the ground state of 
the Hubbard model. The large clusters needed for clar- 
ifying this issue within exact diagonalization methods 
presently render impossible a solution to the problem. 
To estimate the quality of the GA with respect to more 
advanced techniques (which are however restricted to 
smaller clusters and specific boundaries) we have com- 
pared the charge distribution with density matrix renor- 
malization group (DMRG) calculations^ on a 7 x 6 Hub- 
bard cluster doped with 4 holes and for U/t = 12. This 
study also finds metallic stripes. Remarkably the GA 
charge distribution (omitting RPA fluctuations) is almost 
identical to the DMRG. In this case the RPA feedback 
effect is expected to be negligible in the insulator and 
small in the stripe phase since the lowest energy exci- 
tations are gapped in the former and limited to a small 
phase space in the latter. 

The similarity between the exact and the GA charge 
distribution gives further support to the quality of our 
starting point solutions in contrast to HF theory where 
one finds a qualitatively wrong ground state for realistic 
parameters^ In this regard GA+RPA becomes essential 
to proceed. 

The local part of the effective interaction Eq. (0J (low- 
est panel in Fig.^J closely follows the charge modulation. 
This is in contrast to HF where the effective interaction 
between spin fluctuations is uniform and given by the 
bare value of — U. We find that the GA modulated inter- 
action is the dominant effect for obtaining a high-energy 
response in much closer agreement to experimental data, 
as described below, than the related HF+RPA computa- 
tions. 42ii2i2ii The interaction is also strongly renormal- 
ized in our approach, one of the reasons why GA+RPA 
performs much better than HF+RPA upon comparing 
with exact diagonalization results*^ It is interesting to 
notice that the effective value of the interaction we find 
is close to the value usually taken for the bare interaction 
in HF+RPA computations (U w At) in order to obtain 
physically reasonable results4244i2£i This phenomenolog- 
ical decrease of U may compensate for some of the inad- 
equacies of HF+RPA producing a theory more similar to 
GA+RPA. 



IV. RESULTS 

Figures 15161 and 171 report magnetic excitations for var- 
ious dopings and stripe separations. The spectra are 
for selected cuts parallel (right panels) and perpendic- 
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FIG. 5: (Color online) Dispersion of magnetic excitations 
along the cut (jr, 0) — * (-K, 7r) and — > (0,7r) for BC 
stripes and dopings < 1/8 and stripe separation d > 4. 

The dispersion is obtained from 5 contours of uix q (w) (which 
is dimensionless) at regular intervals between and 7. 



ular (left panels) to the stripes and correspond to inten- 
sity plots of the imaginary part of the transverse mag- 
netic susceptibility times the frequency ujx q ( w )- The fre- 
quency factor cancels al/u intensity divergence at the 
Goldstone mode and allows to visualize all energies with 
the same intensity scale. 

The magnetic excitations for SC and BC stripes are 
almost identical with only minor differences regarding 
the intensity distribution and the dispersion as shown 
in Fig. El for rih = 0.13. In the other panels we there- 
fore only report BC solutions which are the ones more 
likely to be stable at low doping as discussed above. This 
similarity is in strong contrast to LSWT calculations for 
which the spectrum is substantially different Indeed, 
within LSWT the number of bands is d for BC stripes 
and d — 1 for SC stripes. In the latter case the core spins 
are assumed to be in S = states and thus to be pas- 
sive. In GA+RPA instead all electrons are active and the 
number of low-energy magnetic bands is d irrespective of 
the symmetry of the stripe. (This is not obvious from 
the figures for large d because the GA+RPA bands have 
small gaps and are only visible in restricted parts of the 
Brillouin zone.) 

As expected a Goldstone mode emerges from the in- 
commensurate wave vector. In GA+RPA a small gap 
(~ 6meV) appears in some runs due to convergence prob- 
lems of the solutions. In these cases we subtract the gap 
to obtain the physically correct gapless spectra. 

The Goldstone mode excitations correspond to long 
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FIG. 6: (Color online) Dispersion of magnetic excitations 
along the cut (tt, 0) — > (vr,7r) and (7r,7r) — > (0, 7r) for BC and 
SC stripes and «/, = 0.13 and stripe separation is d = 4. The 
SC data has been shifted by lOmeV for clarity. We also show 
the fit of the acoustic branch with Eq. Ill 2> . The dispersion 
is obtained from 5 contours of u>x q (w) at regular intervals 
between between and 7. 
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FIG. 7: (Color online) Dispersion of magnetic excitations 
along the cut (it, 0) — > (7r, tt) and (7T, tt) — > (0, 7r) for BC 
stripes at d = 4 and different dopings. The dispersion is ob- 
tained from 5 contours of uJXq ( w ) a ^ regular intervals between 
and 7. 



FIG. 8: Dispersion of magnetic excitations along the cut 
(tt, 0) — > (7r, 7r) and (n,n) — > (0, 7r) for BC stripes oriented 
along the y-direction within LSWT for d — 8 (top) and d = 4 
(bottom). The radius of the points are proportional to the 
intensity times the energy of the excitation. Modes with neg- 
ligibly weight are plotted with a small dot. We used a uniform 
value of the AFM interaction J across all the bonds except for 
the ferromagnetic bonds where the sing is opposite to ensure 
stability and the modulus is Jf = 0.2 J as shown in Fig. 



wave-length transverse fluctuations of the magnetic order 
parameter, i.e. acoustic spin waves. The dispersion rela- 
tion consists of a cone around the incommensurate wave 
vector with an anisotropic intensity distribution. The in- 
tersection of the spin- wave cone with the (q x , u>) -plane 
with q y — tt gives rise to two acoustic branches emerg- 
ing from the incommensurate wave vector Q s . The large 
intensity one disperses towards Qaf = (tTj 7t) where it 
reaches a saddle point (n/i < 0.09) or a local maximum 
(rih > 0.09). We associate the energy at the AFM wave 
vector, ojq, with the resonance energy^ 

The branch emerging from Q s in the cuts and dis- 
persing outwards (i.e. away from Qaf) rapidly looses 
intensity with increasing energy. As can be seen from 
Fig. El this effect is more pronounced for large stripe sep- 
arations. Indeed for = 0.06 the outwards dispersing 
branch is hardly visible. 

In Fig. |5] we show the LSWT results for d — 8 and 



9 



d = 4 BC stripes with Jp = 0.2. We can reproduce 
some qualitative features of the GA+RPA dispersion re- 
lation by taking such a small value of Jf as discussed 
in detail below. The intensity difference effect between 
inwards and outwards dispersing Goldstone mode is also 
observed in the LSWT calculations (Fig. [HJ) , however, the 
difference remains constant as a function of d and is less 
apparent for large stripe separations— i 

Since the outwards dispersing branch has not been de- 
tected yet we cannot make a detailed comparison of the 
relative intensities. Our GA+RPA results suggest that 
this branch should be difficult to observe in the under- 
doped regime. For example for nh = 1/8 the intrinsic 
momentum width of the inward branch observed in the 
experiment will mask the small low energy region where 
the outward branch has significant weight (cf. Fig. 111(1 . 

For rih — 0.16 we expect this problem to be less se- 
vere and the branch should become observable, however, 
low energy INS datai^ do not reveal any signature of the 
outwards dispersing branch. It may be that our theory 
overestimates the corresponding intensity in this case or 
that damping effects, not taken into account in our ap- 
proach, become dominant. More experimental and theo- 
retical work is needed to clarify this point. It should be 
mentioned in this regard that HF+RPA calculations 40 
suggest a further suppression of the intensity of the out- 
wards dispersing modes from the inclusion of a d-wave 
order parameter. 

The two GA+RPA acoustic branches along the 
(ix , % = 1") direction can be fitted by: 

cj(q x ,Tt) =oj \cos(q x d/2)\ (12) 

as shown in Fig. H3 The corresponding velocity perpen- 
dicular to the stripe is: 



The spin wave cone is only weakly anisotropic (cf. 
Figs. Q21 E3J so that the spin-wave velocity along the 
stripe is not far from this estimate. 

In Fig. El we show available experimental data for 
LBCOii n h = 0.125 and LSCO^ n h = 0.16 together 
with the GA+RPA dispersion. There is excellent agree- 
ment between theory and experiment for the position of 
the low-energy acoustic branches and the velocity c of the 
collective mode. We caution, however, that the errors in 
the determination of the experimental spin- wave velocity 
are rather large as shown by the error bars in the right 
panel of Fig. |5J Even a spin-wave velocity independent 
of doping is compatible with the present experimental 
data. To determine the theoretical spin-wave velocity, 
shown also in the right panel, we used Eqs. 1(1211 . 1(13(1 . 
For some runs it was necessary to include a correction to 
Eq. ((12() of the form ui\ sin(q x d)\ with wi in the range 
±6meV. This correction, however, was too noisy due to 
convergence problems of the runs so we include this effect 
as error bars in the theoretical data. 
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FIG. 9: (Color online) Left panel: Dispersion of magnetic 
excitations along the cut (7r,7r) — > (n/4,n) for BC stripes for 
dopings rih ~ 0.1 (short dashed line), 0.125 (full line) and 0.16 
(lon g d ashed line). We include the experimental data from 
Ref. I 111 for rih = 0.125 (triangles) and from Ref. for = 
0.16 (squares). Right panel: velocity of the collective mode 
c divided the lattice spacing a. We report the theoretical 
result (squares) and the experimental result (triangles) both 
obtained by fitting Eq. 1121 to the data. The diamond at 
rih = is the spin-wave velocity in the insulator. 



The theoretical velocity as a function of doping has a 
pronounced minimum at — 1/8 which seems to be 
supported by the available experimental data with the 
mentioned caution on the errors involved. 

Inspection of Eq. 1(13(1 reveals that at small concentra- 
tions the velocity is mainly determined by d ~ 1 /nh due 
to the weak doping dependence of ujq in this regime [cf. 
Fig. If 4( a)], However, for nh > 1/8 the stripe periodic- 
ity stays constant (d = 4) and the rising of qjq (cf. Fig. 
114b ) causes the increase of the velocity and thus the ap- 
pearance of the minimum at nh = 1/8. For nh = 0.1 
the available data for the collective modei^ (not shown) 
imply a velocity substantially larger than what we found 
(as shown in the right panel of Fig. EJ). The predicted 
dispersion is shown in the left panel of Fig. [21 with the 
dashed line. We expect that the saddle point energy uiq 
is similar or slightly smaller than at nh = 0.125 whereas 
the present low-energy experimental data— i extrapolate 
to a significantly larger energy. A high-energy experi- 
ment should clarify the situation on this point. 

For nh > 0.09 the GA+RPA computation additionally 
shows a "roton like" minimum in the direction of the 
stripe (Figs. EE0I- This effect cannot be explained in 
LSWT without invoking longer range interactions. 

The minimum is located close to (tt, 37r/4) which is 
equal to Q s with x and y directions interchanged. It in- 
dicates a tendency of the stripe to develop a spin-density 
wave with period 8a along the stripe. One can spec- 
ulate that for different parameters the system may de- 
velop such an instability resulting in a ground state with 
a checkerboard pattern as observed on the surface of some 
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FIG. 10: Contour plot of \ q (w) as a function of wave vector 
for u> = 45.6meV and rih = 0.125. The plot range for the axes 
is the same as in Fig. 1(c) of Ref. H Wave vectors are in 
units of 2n/a. 



compounds— 

The maximum at Qaf = (t, t) and the roton min- 
imum give rise to a peak in the momentum integrated 
structure factor which is slightly lower in energy than 
the saddle point "resonance" at the AFM wave vector 
CFig. 1161) . Especially at low doping (i.e. large stripe sep- 
arations) the excitations around the (n, tt) region carry 
significant weight which in the integrated spectrum is 
further enhanced due to the saddle-point structure. 

Because of the roton minimum a constant cut at an 
energy moderately below the resonance shows intensity 
both along and perpendicular to the direction of the 
stripe as shown in Fig. Although our study focuses 
on LCO it is remarkable that a similar effect has been re- 
cently observed in an untwinned sample array of YBCO 
by Hinkov and coworkers^ The intensity difference be- 
tween the two directions in also in qualitative agreement 
with our result. 

In Ref. this quasi two-dimensional distribution of in- 
tensity is interpreted as evidence against a "rigid stripe 
array" which is expected to produce a more anisotropic 
pattern. It is clear that Ref. ^| cannot be considered 
as evidence against an oriented stripe array since these 
authors definitively found an asymmetry between x and 
y directions, i.e. a breaking of C4 lattice symmetry. As 
already outlined in the introduction the latter can be 
viewed as the fingerprint of an oriented stripe ground 
state (in the sense of a 'stripe nematic phase—) but may 
also be due to a structural anisotropy as proposed in 
Ref. Without further quantitative estimates one can- 
not distinguish between both scenarios from Ref. 

Our result shows that long-range order with oriented 



stripes is compatible with a quasi two-dimensional dis- 
tribution of intensity between the saddle-point energy 
and the roton minimum. At lower energies we predict 
a strongly anisotropic intensity as long as the disper- 
sion is not cutoff by the spin gap. Can we agree with 
the conclusions of Hinkov and coworkers? This becomes 
a question of semantics. Our stripes are quite different 
from the "rigid stripe array" assumed in LSWT since our 
RPA computation is performed on top of a much softer 
modulation and also has a much larger phase space for 
fluctuations which substantially change the results. In- 
terpreting "rigid stripe array" in the LSWT sense we can 
agree with the statement Ref. [l9| since, as we have shown, 
LSWT in its simpler form is not able to reproduce this 
result. 

We believe that the fact that our results qualitatively 
reproduce the observed two-dimensional intensity map 
observed by Hinkov et al. is evidence for the existence 
of the roton minimum which was obtained^ before the 
data of Ref. became available. Naturally, a quan- 
titative comparison should include consideration of the 
YBCO bilayer structure which would modify the energy 
of the resonance and may affect intensities. In addition, 
a description which allows for fluctuating stripes from 
the starljS may be more appropriate for YBCO mate- 
rials since it can capture the appearance of a spin gap. 
In this regard some details of our map differ from the 
experimental data (round shape, external branches) per- 
haps because of the above reasons. One should keep in 
mind, however, that the map of Hinkov et aZ— is ob- 
tained on the basis of a simple model fitted to 4 scans in 
each direction of which only 3 intersect the ring. Clearly 
details not considered by their model will not show up. 
This should be com pared with a typical experiment at 
ISIS like Ref. UlEci where the maps involve hundreds of 
scans per direction. 

In order to compare with experimental data in twinned 
samples or for systems in which stripes are stacked or- 
thogonal to each other from one plane to the next, one 
should average over the two possible orientations of the 
stripes (vertical and horizontal). Whereas the high- 
energy magnetic excitations in the direction of the stripes 
evolve from the saddle-point at Qaf the spectrum per- 
pendicular to the stripes is composed of several optical 
branches which are not directly connected to the 'reso- 
nance'. Nevertheless it is interesting to observe that in 
GA+RPA the high-energy dispersions parallel and per- 
pendicular to the stripes (Fig. Ill|) are rather close in- 
dicating an almost isotropic magnetic response at large 
frequencies. This high energy 'isotropy' is only present 
up to nh ~ 0.15 where due to the interaction between 
magnetic domains the optical branches perpendicular to 
the stripes start to develop a significant gap structure. 
Fie. II II reveals that also for the optical branches our cal- 
culation yields reasonable agreement with the available 
data— In striking contrast LSWT (Fig. [SJ predicts a 
series of optical branches in the direction perpendicular 
to the stripe that are not observed and an anisotropic 
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FIG. 11: (Color online) Dispersion of magnetic excitations 
along the cut (n, 0) — * (it, n) and (ir, n) — > (0, n) for BC 
stripes at d — 4 and — 0.125. We have averaged the 
data over the two possible orientations of the stripe, vertical 
and horizontal. The dispersion is obtained from 5 contours 
of u>Xq M a t regular intervals between and 10. We also 
show the experimental data from Ref.^3 (courtesy of J. Tran- 
quada). For the latter horizontal "error" bars indicates the 
half width of the peaks whereas vertical "error" bars indicate 
an energy interval of data integration. 



high energy response. The GA+RPA high energy exci- 
tations shown in Fig. 1111 are slightly displaced towards 
higher momenta. However, we find excellent agreement 
for the slope which defines an effective high-energy spin- 
wave velocity. The shift of our excitations towards higher 
momenta may be due to the use of the oversimplified one- 
band Hubbard model. It would be interesting to see how 
the high energy part of the dispersion appears in the 
three-band Hubbard model. 

A quantitative fit to the X-type dispersion shown in 
Fig. ^]can also be achieved with coupled spin-leg ladder 
models by fine tuning the interstripe coupling to a quan- 
tum critical point.— These theories, however, rely on an 
even charge periodicity of the stripes (as for rih = 1/8). 
In contrast, we obtain qualitatively similar spectra for 
nh = 1/8 (d = 4) and rih = 1/10 (d = 5) in accord with 
experiment in LCQiiii On the other hand our approach 
cannot be easily generalized to incorporate spin gapped 
states as in YBCO which are naturally explained in a 
ladder picture^^LSi^ Another important difference of 
the ladder scenario with respect to the present work is 
that at high energies the ladder computation finds prop- 
agating modes only in the direction of the stripe whereas 
our calculation yields more isotropic spectra in agreement 
with recent experiments^ 



Constant energy cuts of the dispersion change shape 
as the energy of the cuts is increased. We have already 
shown 1 that our data closely resemble the changes seen 
in the experiment at rih = 1/8. 11 In Figs. 1121131 we 
show contour plots of the transverse susceptibility for 
a strongly underdoped (rih = 0.08) and an almost opti- 
mally doped (n/i = 0.16) system averaged over the two 
possible orientations of the stripes. In the first case we 
consider an underlying d = 6 BC stripe solution which 
yields a saddle-point energy at ujo = 43meV. Below this 
energy the 'resonance feature' splits into the four incom- 
mensurate Goldstone modes (from the superposition of 
vertical and horizontal stripes). The associated spin- 
wave cones display a strongly asymmetric intensity distri- 
bution where the weight is mainly confined in the parts 
closest to Qaf in accord with the intensity difference 
of the acoustic branches shown in Fig. [5] and discussed 
above. 

Above too the intensity spreads into a ring shaped fea- 
ture which increases with energy. The strongest scatter- 
ing is on the corners of a square at 15meV and rotates by 
45° at 60meV resembling the behavior found at rih = 1/8 
in theorjii and experiment. 

The case of d = 4 BC stripes at doping rih = 0.16 
is shown in Fig. 1131 The saddle-point energy is at 
u)q = 102mcV. At low energies the spin-wave cone dis- 
plays significant weight in all directions which leads to the 
pronounced intensity of the outwards dispersing branch 
as discussed already in context of Fig. [7| The complex 
structure at 87meV is composed of the incommensurate 
Goldstone modes and the branch dispersing from Qaf 
towards the roton minimum which constitute the smaller 
ring around Qaf- Additionally one can observe a larger 
ring which intensity comes from the branch dispersing 
away from the roton minimum along the stripe direction. 
A weaker feature along the diagonals originates from the 
outwards dispersing branches. Above the saddle-point 
energy the intensity distribution again forms a ring which 
increases in diameter with increasing energy. 

Fig. [21 displays the doping dependence of the reso- 
nance frequency and intensity for both BC and SC stripes 
evaluated within GA+RPA. In the strongly underdoped 
regime (0.05 < rih < 1/8) we find an approximately lin- 
ear relation fkoo rs a + firih with (3 w 250meV and a w 20 
(30) meV for BC (SC) stripes in the GA+RPA. Note that 
for infinitely distant stripes one expects u>o — > so that 
the offset energy a is an artefact of the linear approxi- 
mation in the considered doping range. A qualitatively 
similar linear behavior has been found in Ref. lo2t 

We also show in Fig. I14f a1 (full squares) the dop- 
ing dependence of the resonance energy evaluated within 
LSWT. We assumed that the doping is given by rih — 
l/(2d) which is valid for half-filled stripes with rih < 1/8. 
The energy of the resonance scales in a similar way in 
both GA+RPA and LSWT for n h < 1/8 despite the fact 
that we take the values of J and Jp as doping indepen- 
dent. This is not surprising since for rih < 1/8 the charge 
profile of the stripes do not overlap and doping proceeds 
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FIG. 12: Contour plots of Xq ( w ) as a function of wave vector 
at selected energies for doping nn = 0.08. The underlying 
GA solution is a d — 6 BC stripe array. The wave vector 
coordinate system has been rotated by 45° in such a way that 
the magnetic Brillouin zone of an antiferromagnet occupies 
all the window of the plot with Qaf ~ (tt, tt) at the center. 
Wave vectors are in units of 2ty/(^/(2)o). 



by varying the distance among the stripes keeping the 
electronic structure on the core of the stripe essentially 
unchanged Therefore in this respect it makes little 
difference at low energy to consider the stripes as local- 
ized objects as in LSWT with a fixed value of Jf or as 
more extended objects as in GA+RPA. It is interesting 
that the effective value of Jf in order to have an over- 
all agreement with the GA+RPA results and with the 
experiment is quite small as found in other LSWT stud- 

41 

lesiai 

For rih > 1/8 effects due to the overlap of the stripe 
cores become evident. In our previous computationsS^i 
within the GA we have found that in this regime it is 
energetically favorable to add holes into the d = 4 stripe 
structure consistent with the doping independent incom- 
mensurability seen in INS experiments^ For both BC 
and SC stripes this leads to a strong increase of the reso- 
nance energy and therefore to a deviation from the linear 
ujq vs. nh behavior at low doping. An alternative sce- 
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FIG. 13: Contour plots of Xq ( w ) as function of wave vector at 
selected energies for doping rej, = 0.16. The underlying GA 
solution is a d = 6 BC stripe array. The wave vectors units 
are as in Fig. 1121 



nario is phase separation as explained below. 

The resonance mode is strongly influenced by the ef- 
fective magnetic interactions across the domain wall. For 
example, in case of a SC stripe depicted in Fig. 1151 the 
spins adjacent to the core site are bridged by an effective 
AFM interaction J' as it is assumed in LSWT models. 
The energetics of the resonance mode can be easily under- 
stood by adding a small transverse component with mo- 
mentum Q = (tt, tt) to the domain wall structure which 
provides a simplified picture for the quasi-classical dis- 
placements of the spin in the resonance mode as shown 
in Fig. 1151T b') . The spins in the AFM regions keep their 
relative angle constant whereas the angle formed by the 
spins across the domain wall strongly fluctuates, giving 
the dominant contribution to the energy of order J' . A 
similar picture holds for BC stripes with Jp at the place 
of J'. (In reality the resonance mode is not completely lo- 
calized in the domain wall, as suggested by Fig. ED which 
makes the dependence on J' of Jf more complicated.) 

The strong increase of the resonance energy with dop- 
ing for rih > 1/8 [Fig. I14f a')] can be understood if the ef- 
fective value of Jf (J') also increases across the BC (SC) 
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FIG. 14: (color online) Doping dependence of the resonance 
frequency and intensity for BC (left panels) and SC (right 
panels) stripes within the GA+RPA. The circles for n/, > 1/8 
correspond to d = 4 solutions whereas the triangle corre- 
sponds to the d — 3 solution. The full diamond symbol is 
the experimental result from Ref. ITU The full squares in (a) 
are the LSWT result using J F = 0.2J and J = lOOmeV. 



domain walls. As shown in the Appendix B these couplings 
are dominated by transverse hopping processes of doped 
holes. The increase of the stripe filling from v = 1/2 for 
nh > 1/8 induces a concomitant enhancement of Jf {J') 
explaining the behavior of the upper panels of Fig. 1141 

Since upon doping the saddle point energy for nh < 
1/8 increases more slowly than the experimental T c our 
calculations predict a non-constant ratio between loq and 
kT c in the underdoped regime of LCO. This ratio de- 
creases from uJo/kT c w 27 for rih = 0.05 to Lu /kT c w 17 
for rih — 1/8. This is significantly larger than the reso- 
nance energy to kT c ratio ~ 5.4 as obtained for YBCO. 
Therefore if the saddle point excitation is important for 
the superconducting mechanism its effect on T c is mod- 
ulated by other factors. 

In LCO materials the only experimental value for loq 
available up to nowii (shown as the full diamond in 
Fig. ll4f> falls exactly on the curve for BC stripes whereas 
it is slightly below the Uq(1/8) result for SC ones. Al- 
though the difference is too small to be conclusive (given 
the approximations involved) the fact that BC stripes 
give the best fit is in agreement with our previous find- 
ing within the three-band model^ that BC stripes are 
the more stable textures in the underdoped regime of 
LCO cuprates. 

One can obtain two more points for the resonance en- 
ergy by extrapolating the experimental data of Ref. ^1 
with the aid of Eq. I|12|l . One obtains a similar result as in 
the right panel of Fig. 03 The extrapolated experimental 
resonance is in good agreement with theory for nh = 0.16 
whereas it is at much higher energies for = 0.1. 

Around optimal doping one finds a proliferation of 
mean-field solutions with energy similar to stripes prob- 
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FIG. 15: a) A schematic cut perpendicular to the stripe 
through the Neel-type spin structure of d — 4 SC stripes, b) 
Resulting spin structure upon adding transverse spin compo- 
nents S± ~ cos(QR) with Q = (ir, tt) 



ably indicating a melting of these textures. Therefore in 
this regime the nature of the ground state is less clear. 
Restricting to vertical stripes we find a transition towards 
d = 3 stripes both in the one— and three-bandSi model 
at rih — 0.19 ~ 0.24. Experimentally the incommensu- 
rability e for doping > 1/8 has been determined for 
LSCO 10 and LCO codoped with Nd 9 . In the former ma- 
terial e s» 0.125 up to nh = 0.25 whereas for the latter 
compound the incommensurability reaches e ~ 0.15 for 
nh = 0.2 which is compatible with d = 3 stripes. For this 
solution the resonance frequency and intensity is shown 
with a triangle in Fig. 1141 Clearly more measurements of 
the dispersion and in particular the resonance at doping 
n/j > 1/8 will clarify the situation. 

In Fig. E| we show the density of magnetic excitations 



(14) 



Here the sum is restricted to the magnetic Brillouin zone 
for commensurate AF order to conform to the usual ex- 
perimental practice and N denotes the number of lattice 
sites. 

As doping increases spectral weight is transfered from 
high energy to the resonance region in good agreement 
with experiment ifiLS The increase of weight around the 
resonance energy is not a matrix element but a density 
of states effect. Indeed the intensity at the resonance 
decreases with doping (Fig. 114(1 which, however, is over- 
compensated by the strong Van Hove singularity pro- 
duced by the saddle-point like dispersion around (7r,7r). 
We have performed preliminary calculations of S(u>) for 
nh < 0.05, where stripes are diagonal, which in fact in- 
dicate that at very low doping the density of magnetic 
excitations acquires again the characteristics of the un- 
doped AFM. 

For rih > 0.09, both the roton minimum and the 
saddle-point ('resonance') contribute to the peak. The 
corresponding energies are marked as vertical bars in the 
panels of Fig. ^] At nh = 0.1 both energies almost co- 
incide and give rise to the peak in S(u>) at w = 45meV 
whereas for higher dopings they are clearly separated. 
At nh — 0.125 the overall structure of S(ui) consists 
of a dominant peak around the 'resonance' frequency 
and several humps at higher energy arising from the op- 
tical bands. We obtain good agreement with the ab- 
solute intensity of the experiment Ref. by choosing 
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FIG. 16: Density of magnetic excitations S(ui) (see text) for 
three doping levels n h = 0.083,0.1,0.125,0.132. Theoretical 
curves include a factor ZdT]± = 0.25. The vertical bars indi- 
cate the energies of the saddle-point at Qaf and the roton- 
type minimum in the dispersion along the stripes. Experi- 
mental data for nn = 0.125 are from Ref. I 111 

ZdV± ~ 0-25 (cf. Eq. 1111) ) which is in good agreement 
with the estimates of Ref. |6jJ 



V. CONCLUSIONS 

In this paper we have studied the doping dependence 
of spin excitations in high-T c materials of the LCO type 
based on a stripe picture. The necessity for an under- 
standing of these excitations is suggested by the direct 
relation found between incommensurate scattering par- 
allel to the Cu-0 bond and superconductivity in LCO 
compounds. 14 i 15 

The comparison with the limited experimental data 
available so far gives good agreement for both acoustic 
low energy magnons and the high energy optical spin 
excitations. Also the absolute intensity of the modes is 



in reasonable agreement with experiment taking into ac- 
count the results of Ref. |(37]- The fact that we can rec- 
oncile our computations with experiments is remarkable 
insofar given the simplicity of the hamiltonian and the 
fact that the three model parameters U, t and t' are 
fixed by strong constraints previous to the computations 
of spin excitations in the doped phase. 

We provide a non trivial prediction for the dependence 
of the velocity of the acoustic spin-wave mode and the 
energy of the saddle point on doping. The spin-wave ve- 
locity has a dip at rih — 1/8 which coincides with a sim- 
ilar anomaly seen in the doping dependence of T C M&^ 
Although the superconducting anomaly is probably re- 
lated to charge ordering along the stripe, the coincidence 
is intriguing. 

The nature of the ground state is not so clear as one 
approaches the overdoped regime since one has to take 
into account the melting of the stripes and non-Gaussian 
fluctuation effects beyond RPA will become important. 
Also it is conceivable that phase separation between the 
underdoped stripes and an overdoped Fermi liquid occurs 
as mentioned in Ref. l50l Since we assume a uniform 
stripe phase one can find evidence in favor or against 
this scenario by studying the doping dependence of the 
saddle point energy. In the phase separation scenario we 
would expect that at some point the population of the 
domain walls by holes saturates thus hampering a further 
increase of the saddle point energy with doping. 

The GA+RPA results considerably differ from the pre- 
dictions of LSWT, however, it is interesting that some 
features can be qualitatively reproduced by this sim- 
ple theory. These include: i) the difference in intensity 
among the two acoustic branches emerging from the in- 
commensurate point, 81 ii) the strong intensity at the sad- 
dle point associated with the resonance, Hi) the fact that 
the first optical branch in the direction perpendicular to 
the stripe is separated from the acoustic branch at (tt, 7t), 
and iv) the doping dependence of the resonance energy 
for rih < 1/8 as shown in Fig. 1141 

Regarding the differences between LSWT and 
GA+RPA the most notable ones are the absence of the 
roton minimum in LSWT and the different number of 
bands among SC and BC dispersions in LSWT vs. the 
similarity of the dispersion for the two symmetries in 
GA+RPA. Moreover we have found that the dispersion 
of the optical branches perpendicular to the stripes sig- 
nificantly differs between GA+RPA and LSWT. For ex- 
ample, upon comparing Figs.[7]and|Hlit turns out that the 
first optical branch is inverted with respect to its mean 
energy. 

Within LSWT the coupling between adjacent AFM do- 
mains is described by an effective ferromagnetic exchange 
Jp in case of BC stripes (or an AFM coupling J' for 
SC textures). It is interesting that the effective value of 
Jf (J'), in order to have an overall agreement with the 
GA+RPA results and with experiment, is quite small as 
found in other LSWT studies. 41 With this small value one 
may wonder whether the magnetic system has long-range 
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order at all or whether it is disordered due to quantum 
fluctuations^ In any case we expect the system to be 
close to a magnetic quantum critical point. This is in 
agreement with the experimental results of Aeppli and 
collaborators which indeed see nearly quantum critical 
behavior close to the incommensurate wave vectors 

Regarding the RPA approach that assumes a uniform 
Fermi liquid as a starting point 27 i 28 i 29 i 30 i 31 i 32 i 33 i 34 i 35 it 
has the advantage that it does not break C4 lattice sym- 
metry nor SU(2) spin rotational symmetry in agreement 
with what one expects in overdoped materials. On the 
other hand since the starting point is a Fermi liquid this 
approach can not be pushed to the dilute regime or to en- 
ergies of the order of the charge transfer gap, i.e. describe 
charge excitations across the Hubbard bands. Within 
this scheme one obtains a dispersion relation for spin ex- 
citations that differs considerably from ours. For exam- 
ple, there is no optical branch that continuously merges 
with the saddle pointi 3 ^* Also the low energy response 
never separates into four spots in momentum space as the 
spin-wave cones for a striped system. Instead the incom- 
mensurate spin response forms a two-dimensional struc- 
ture at very low energies. Also the integrated spectral 
weight is concentrated in a much smaller energy ranged 
Finally for untwinned systems spectral weight is expected 
to appear in both directions with some modulation due 
to structural anisotropy 3 - whereas in our theory at the 
lowest energies practically all the weight is in the direc- 
tion perpendicular to the stripes. 

The natural question which arises is whether the re- 
sults of the present investigation may also be relevant 
for an understanding of spin excitations in YBCO. The 
acoustic dispersion shows a low-energy spin gap (A ~ 
30mcV at optimal doping) and thus does not reach the 
incommensurate wave-vector at u> — 0. In this case it is 
more likely that the system is in a quantum disordered 
spin phase as suggested by the ladder theoriesi 36 i 37 i 39 
One expects that the systems shows short range order 
with a correlation length of the order of A /(he) ~ 5a 
and that for energies larger than A the system resem- 
bles that of an ordered phase. For this reason one can 
expect that a computation as the present one is suit- 
able for a description of the universal high energy spin 
response ^LiiSS One encounters a similar situation in the 
one-dimensional Hubbard model for large U where a com- 
putation of RPA fluctuations on top of a spin-density 
wave will in many respect be in better agreement with 
the exact solution (which has no long-range order) than 
RPA fluctuations on top of a homogeneous Slater deter- 
minant. Alternatively a scenario of fluctuating stripes 
where also the charge loses its long-range order can cap- 
ture the effect of a spin-gap^ One has to keep in mind, 
however, that too much orientational fluctuations can re- 
sult in an isotropic state which is not what is found in 
the experiment of Ref. 

In YBCO one finds an apparent correlation between 
the resonance frequency and T c in the sense that both 
display a similar parabolic shape as a function of doping 



so that loq w 5.3kT c . Measurements of luq have been per- 
formed from 7% underdoping up to 2% overdoping (see 
e.g. Refs. HH). Due to the spin gap it is more difficult to 
extract the behavior of the incommensurability from the 
data. The authors in Ref. 16 find that e already saturates 
for doping w 0.1 compatible with a stripe spacing of 
d = 5 (extrapolation to uj — implies a larger incom- 
mensurability similar to LCO). It is then reasonable to 
assume that the saturation corresponds to the situation 
where holes are doped into a stripe structure with fixed 
distance as in case of LCO for rih > 1/8. 

We find that the resonance energy as a function of dop- 
ing tends to have a saturation or a slower rate of growth 
with doping in the underdoped and the overdoped regime 
[Fig. H3k]. This is compatible with Fig. 29 in Ref. El 
which reports a flattening of the luq vs. T c relation for 
both optimally doped and strongly underdoped samples. 
In the latter case this would correspond to the transition 
towards the regime where stripes are created at distances 
d>5. 

Finally our investigations suggest an explanation for 
the apparent two-dimensional magnetic scattering be- 
low the resonance in detwinned samples of YBCO^i 9 " Ar- 
guments against a "rigid" stripe in this compound are 
based on the assumption that in a detwinned material 
stripes preferably align in the direction of the CuO chains 
running along the b* direction. From spin-only stripe 
models 36,37,41 '— one then expects a low energy magnetic 
response only in the a* direction. Instead, the experi- 
mental datai&i 9 " show an anisotropic intensity modula- 
tion with maxima along the a* direction but also signifi- 
cant weight along b* . Lateral fluctuations of the stripes^ 9 , 
can also explain this result which in our harmonic RPA 
approach are presenljS but decoupled from the magnetic 
modes. 

On the other hand our present computations show 
that the observed two-dimensional, but anisotropic in- 
tensity modulation can be explained by the roton like 
minimum in the dispersion. In addition our spectra dis- 
play an almost isotropic response for energies well above 
the saddle-point in agreement with neutron scattering re- 
sults on a detwinned and underdoped YBCO sample by 
Stock and coworkers 22> The calculation of magnetic exci- 
tations from stripes within the time-dependent GA there- 
fore yields a non-trivial crossover from one-dimensional 
behavior at low energies to a quasi two-dimensional re- 
sponse at high energies with an intermediate anisotropic 
but two-dimensional structure at the roton minimum. In 
contrast spin-only models yield a one-dimensional spec- 
trum over the whole range of energies. 

After this work has been submitted and posted^i new 
INS data for an untwinned sample array of YBCO be- 
came available^ 3 - In our opinion the strong anisotropy 
found at the lowest energies is incompatible with uni- 
form Fermi liquid theories. In contrast some of our pre- 
dictions based on stripes are qualitatively reproduced. 
Above the resonance the intensity is similar for both di- 
rections in qualitative agreement with Figs. 1516*171 and in 
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disagreement with the ladder theories ^L2L22i Moderately 
below the resonance there is intensity in both directions 
whereas at lower energies the intensity becomes much 
more anisotropic [c.f. Figs. 2(a) and (b) of Ref. 93] which 
is qualitatively consistent with a roton minimum. On 
the other hand our computations do not include super- 
conductivity, temperature effects and damping which are 
clearly important in the experiments. In addition specific 
computations for bilayer systems are in process in order 
to investigate the relevance of our results for YBCO ma- 
terials more quantitatively. 
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APPENDIX: DOPING DEPENDENCE OF THE 
MAGNETIC COUPLING ACROSS STRIPES 

In this Appendix we present simple arguments to show 
how the filling of the domain walls by holes for nh > 
1/8 affects the effective magnetic interaction across the 
stripe. Our arguments will be based on strong-coupling 
and we decouple the motion of charge carriers perpen- 
dicular and along the domain walls. The kinetic energy 
gain of the latter is responsible for the stability of v = 1/2 
half-filled stripes. 51 In the following v is fixed and we con- 
centrate on the transverse motion which gives rise to the 
magnetic coupling between antiphase AFM domains. 

SC stripes: For half- filled stripes (n^ < 1/8) the core 
will be essentially in a mixture of configurations of the 
kind: 

T i (50%) 
T T I (25%) 
T I I (25%) 



where we have represented three sites perpendicular to 
the stripe the central one corresponding to the core. 
denotes a hole and we focused on three sites in which 
the prevailing configuration of the external spins is f and 
|. Since for half- filled stripes a hole occupies a core site 
with probability 50%, the absence of a hole implies the 
presence of an up or down spin with probability 25% re- 
spectively. One can draw a similar picture in the case in 
which the prevailing configuration of the external spins 
is i and f. The hole in the core favors the AFM cou- 
pling of the adjacent spins through virtual processes of 
the kind jj where jj denotes a doubly occupied site. 
Therefore the gradual population of the core with holes 
for rih > 1/8 should enhances the coupling J' thus lead- 
ing to an increase of the resonance energy as discussed in 
Sec. HV| 



BC stripes: Stability of the structure requires that the 
core is composed of two sites which are coupled by an 
effective ferromagnetic interaction Jp. The core spins 
are coupled by an effective antiferromagnetic exchange J 
to the external spins. For the case in which the prevailing 
configuration of the external spins is j and j (like in the 
upper row of Fig. ^) one has for a half-filled stripe the 
typical configurations: 



I T I (25%) 
4 T I (25%) 

I T T I (50%) 



The first two configurations favor the antiphase align- 
ment through processes of the kind || J. and 
J. T-l 0. Since these process involve holes one ex- 
pects an increase of the effective Jp when the numbers of 
holes in the domain wall starts to deviate from half-filling 
for n h > 1/8. 
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